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CHAPTER    1 


INTRODUCTION 


The  radio  index  of  refraction  in  the  atmosphere 
is  defined  as  the  ratio  of  the  speed  of  light  to  the  velocity 
of  a  radio  wave  in  the  atmosphere.   There  are  two  methods 
of  measuring  the  index  of  refraction,  the  direct  method  em- 
ploying radio  ref ractometer s  and  the  indirect  method  utiliz- 
ing standard  weather  observations.   The  indirect  method  re- 
quires the  measurement  of  temperature,  pressure,  and  humidity 
with  the  subsequent  conversion  to  the  index  of  refraction 
using  the  formula  below: 


(n-1)  X  106  =  77-6P  "  lie  +  5.8  X  105 

T        T  T* 


(Bean  et  al,  1966) 


or  a  more  convenient  expression  in  N  units, 


N  =  (n-1)  x  106  =  77-6P  "  lie.  +  5-8  ><  1°5 

T        T         TH 


where,  n  =  the  index  of  refraction  of  the  atmosphere 


p  =  barometric  pressure  in  millibars 

T  =  temperature  in  degrees  Kelvin 

e  =  water  vapor  pressure  in  millibars 

The  operation  of  radio  ref ractometer s ,  used  for 
direct  measurement,  is  based  on  the  fact  that  the  resonant 
frequency  of  a  cavity  resonator  is  inversely  proportional 
to  the  index  of  refraction  of  the  gas  filling  the  cavity. 
Most  ref ractometer s  utilize  two  resonant  cavities,  a  sealed 
reference  cavity  and  a  ventilated  sampling  cavity.   Since 
the  reference  cavity  is  sealed  and  maintains  a  constant  res> 
onant  frequency  (within  the  limits  of  small  temperature 
and  pressure  changes),  the  change  in  resonant  frequency  of 
the  sampling  cavity  due  to  the  change  in  the  index  of  re- 
fraction of  the  air  passing  through  it  can  be  detected. 
There  are  a  number  of  ways  of  detecting  this  resonant  fre- 
quency change  and  converting  it  to  the  index  of  refraction, 
for  example,  the  Crain  type  ref ractometer  measures  the  dif- 
ference in  frequencies  and; 


Af 


=  -An,  AN  =  An  x  10* 


The  Birnbaum  type  ref ractometer  measures  the  phase  differ 
ence  between  the  two  cavities. 


Considering  only  the  accuracy  of  the  measurement 
and  the  time  required  to  evaluate  the  data,  the  direct 
method  is  desirable  since  the  value  of  N  is  determined 
directly  and  only  one  sensor  is  required  instead  of  three. 
However,  the  disadvantages  of  high  cost,  complexity,  and 
maintenance  required  outweigh  the  advantages  and  the  majority 
of  refractive  index  measurements  to  date  have  been  made  in- 
directly. 

The  ideal  situation,  especially  in  the  field  of 
radio  propagation  studies  where  precise  measurements  of  the 
index  of  refraction  are  required,  would  be  to  have  a  simple 
inexpensive  method  of  measuring  the  index  of  refraction  di- 
rectly.  One  possible  method  of  accomplishing  this  is  to  use 
the  backscattering  characteristics  of  a  ventilated  resonant 
cavity  as  the  cavity  falls  through  the  atmosphere.   As  the 
cavity  falls  through  the  atmosphere,  it  can  be  tracked  by 
radar  and  the  change  in  backscattering  due  to  a  change  in 
the  resonant  frequency  can  be  detected  and  related  to  the 
index  of  refraction  of  the  gas  passing  through  the  cavity. 

In  recent  years  a  considerable  amount  of  research 
has  been  directed  to  the  reduction  of  backscattering  or 
radar  cross -section  by  impedance  loading  particularly  in 
the  resonant  region.   Chen  and  Liepa  (1964)  studied  the 


scattering  cross -section  of  a  thin  center -loaded  cylinder 
and  demonstrated  the  reduction  of  the  radar  cross -section 
due  to  loading.   Liepa  and  Senior  (1966)  studied  the  scat- 
tering behavior  of  a  metallic  sphere  with  circumferential 
slot  loading.   Plonus  (1967)  investigated  the  scattering 
behavior  of  a  spherical  shell  with  a  finite  slot  for  loading. 
Chang  and  Senior  (1969)  studied  the  scattering  behavior  of 
a  thin  spherical  shell  with  a  circular  aperture  and  demon- 
strated that  there  is  a  large  reduction  of  the  radar  cross- 
section  in  the  vicinity  of  the  cavity  resonant  frequency. 

This  thesis  is  concerned  with  determining  if  the 
resonant  cavity  approach  to  the  direct  measurement  of  the 
radio  index  of  refraction  in  the  atmosphere  is  feasible. 
A  thin  spherical  copper  shell  with  circular  apertures  has 
been  selected  for  this  study  because  of  the  symmetry  pro- 
perties of  the  sphere  and  the  circle. 


CHAPTER         II 


THEORY 


2.1   Solutions  to  Wave  Equation  in  Spherical 
Coordinates  with  Axial  Symmetry 

With  axial  symmetry  Ramo,  Whinnery  and  Van  Duzer 

(1965)  approaches  the  problem  by  separating  the  solutions 

into  waves  with  components  Er,  Eq,    and  H^  and  waves  with 

components  Hr,  Hq,  and  E<j>.   Using  the  spherical  surface 

(r  =  constant)  as  the  transverse  surface  the  wave  with 

components  Er,  Eg,  and  H$  is  called  the  TM  type  and  the 

one  with  components  Hr,  He,  and  E  <j>  is  called  the  TE  type. 

Considering  only  the  TM  type  with  H   =  0,  H   =  0, 

r  B 

E^    =    0    and    h/h$    =    0,    Maxwell's    equations    in    spherical 
coordinates    become: 

VX    H=ja>€    E 


a  d(H       sin    9)  a 

r  » '   +     e 


r    sin    6 


V    X    E    =    -    joouH   =    a 


4> 


Be 


d(rEft) 


-d(rH    )" 


br 


c)r 


Be 

3 


(1-1) 


(1-2) 


then  since  H  =  a  ^H .  and  E  =  arEr  +  aeEQ 


1    5(rEe)  3Er 

F  — §J F  5e"  =  -jm^H« 


r    sin    © 


cKH^    sin    e) 
d© 


JOJ    €     Er 


.1  ^v 

7  -5 =  ^  eEe 


(1-3) 


(1-4) 


(1-5) 


solving  for  E   and  E„  in  equations  (1-4)  and  (1-5)  we  get 


E   = 


E.  = 


1 


d(h\  sin  9) 


jo)  e  r  sin  ©     d© 

1     d(rH^) 
jco  e  r    cir 


(1-6) 


(1-7) 


substituting  equations  (1-6)  and  (1-7)  into  equation  (1-3) 
and  carrying  out  the  required  operations,  we  get: 


a  (rHJ      1  a 


!    d(H$  sin  ©) 


£r2 


+  —  — - 
r 


sin  ©      d© 


be 


J  +  oo  u  e  rH.,.  =  0      (1-8) 


Letting  k   =  od  ue  (Ramo  et  al.  ,  1965)  solved  the  second 
order  differential  equation  (l-8)  by  the  product  solution 
method  and  found: 


H*  = 


—  J  ,    i   Jkr) 
£T  n+(l/2)' 


B 


^ 


/,  n+(l/2) 


(kr) 


p  (cos  e) 

n 


(1-9) 


where,  Jn+(l/2)  =  a  Bessel  function  of  the  first  kind 
^n+(l/2)  =  a  Bessel  function  of  the  second  kind 
Pn        =  a  Legendre  polynomial 
If  the  problem  contains  the  origin  Bn  must  equal  zero  be- 
cause Nn  +  (W2)  is  infinite  at  the  origin.   If  the  problem 

extends  to  infinity  the  Bessel  functions  must  be  combined 

(2) 
into  the  second  Hankel  function,  Hn+(]/2)  =  Jn+(l/2)  " 

j^n+(l/2)  "k°  represent  the  wave  traveling  away.   Ramo  et  al 
(1965)  uses  the  notation  Zn+(]/2)  to  represent  the  particu- 
lar combination  of  Jn+(i/2)  and  Nn+(l/2)  required  for  the 
problem  under  consideration  and  then: 


A    _ 

HA  =  ~  P  (cos  9)Z   ,„  /ox  (kr) 
nr      nv      '    n+(l/2) 


(1-10) 


Substituting  equation  (1-10)  into  equations  (1-6)  and  (1-7) 
and  differentiating,  we  get: 


A  P1(cos  6) 


E„  = 


n  n 


jo)  €  r 


3/2 


nV(l/2)(kr)  "  krZn-(l/2)(kr) 


(1-11) 


8 


E   = 
r 


A  nZ   ,  .  .  (kr) 
n   n+(l/2)v 


jcn  e  r 


3/2 


sin 


cos  e  pi (cos  e)  -  p^+1(cos  e) 


(1-12) 


The  same  approach  can  be  used  to  determine  the 
field  equations  for  the  TE  mode  where  B/c)<t>  =  0,  Er  =  0, 
Ee  =  0,  and  H<|)  =  0  or  by  the  principal  of  duality  replac- 
ing Er  by  Hr,  Ee  by  Hq  ,     H<j)  by  -  E$,     £  by  \i ,     and  using  a 
new  constant  B   in  equations  (1-10),  (l-ll),  and  (1-12). 
Applying  the  principal  of  duality  the  following  equations 
are  obtained  for  the  fields  of  the  TE  mode: 


B 


E*    -^_P^(cos  e)Zn+(i/2)(kr) 


-B  P1(cos  e) 
HQ  =  _n  n 


jcDiar 


3/2 


nZn+(l/2)(kr)  -  krZn_(l/2)(kr) 


(1-13) 


(1-14) 


B  nZ   ,   .  . (kr) 
Hr  =   "   n+(l/2) 

jcD|j.r3/2  sin  9 


cos  e  p^(cos  e)  -  pJ+1(cos  e) 


(1-15) 


2.2  Spherical  Cavities 


The  TE,  ni  and  TM, fi,  modes  of  a  spherical  cavity 
will  be  analyzed  using  equations  (1-10)  through  (1-15). 
Assume  the  cavity  to  be  closed  with  inner  radius  a  and 


that  the  shell  is  perfectly  conducting.   The  Bessel  func- 
tion of  the  second  kind  must  be  eliminated  from  our  equa- 
tions since  the  origin  is  included  in  the  area  of  interest, 
therefore,  the  field  equations  for  the  TE]_oi  mode  in  the 
cavity  are : 


B 
E<,  =  -kPi(cos  e)j3/2(kr) 
/r 


(2-1) 


B  P1(cos  e) 


He  = 


1  1 


joour 


3/2 


J3y2(kr)  -  krJj/gCkr) 


(2-2) 


Hr  = 


V3/2(kr) 


jcour 


3/2  sin  6  _ 


cos  e  pJ(cos  e)  -  pJ(cos  e) 


(2-3) 


where  pj(cos  9)  =  sin  0 

pl(cos    9)    =    3    sin    9    cos    9 


J3/2(kr)    = 


rtkr 


'sin    kr 


-    cos    kr 


kr 


J\  /g(kr)    =/ sin    kr 


rtkr 


Making  the  substitutions  indicated  above  in 
equations  (2-1)  through  (2-3)  we  get  the  following  equations 
for  the  TE-|_oi  mode  in  the  cavity: 


10 


B   sin  9   2 


E  *  = 


"sin  kr 


kr 


cos  kr 


-B1    sin  0  /  2 


kr  sin  kr  -  cos  kr 


-  cos  kr 


jo>|ir 


(2-4) 


(2-5) 


(2-6) 


Since  we  made  the  assumption  that  the  shell  of  the  cavity 
was  perfectly  conduct ing,  the  tangential  electric  field  at 
r  =  a  must  go  to  zero,  therefore: 


E.  =  0  =  -i 


B_  sin  e  /  2 
ink 


sin  ka 


ka 


-  cos  ka 


0  =  sin  ka  -  cos  ka 
ka 

0  =  tan  ka  -  ka 

tan  ka  =  ka 

Tan  ka  =  ka  is  the  condition  for  resonance  for 
the  TE]_Q2.  m°de  in  a  perfectly  conducting  closed  cavity. 
The  solution  of  this  transcendental  equation  can  be  solved 
graphically  and  the  solution  is  ka~4.5. 

For  the  TM-lq-l  mode  we  have: 


Ha  — 


A   sin  9  /  2 


sin  kr 


-  cos  kr 


kr 


11 
(2-7) 


Eo  = 


-    A       sin    9  J2 
Jcjo    e    r      \l  Jtk 


-  kr  sin  kr  -  cos  kr 


2A   cos  6  |  2 


Er  =- 


joo  e  r   v  ^k  [   kr 


sin  kr 


cos  kr 


(2-8) 


(2-9) 


As  for  the  TE-^q-j_  mode  the  tangential  electric 
field  must  equal  zero  at  the  shell  for  the  resonant  condi- 
tion for  the  TMjloi  mode,  therefore: 


0  =  sin  ka  _1  -  ka  -  1 
cos  ka  ka 


tan  ka  = 


ka 


ka 


ka 


1  -  (ka) 


The  graphical  solution  to  this  equation  is 
ka%2. 75. 

The  figure  of  merit  or  "Q, "  of  the  cavity  is  de- 
fined as  the  ratio  of  the  energy  stored  at  resonance  to 
the  energy  loss.   Ramo  et  al.  (1965)  defines  the  Q,  for  the 
TM-,q-]_  mode  as  approximately  equal  to  the  ratio  of  the  in- 
trinsic impedance  of  the  gas  filling  the  cavity  to  the 


12 

surface  resistivity  of  the  metal  in  the  shell.   Then  for  a 
copper  cavity  filled  with  free  space  the  Q,  would  "be. 


Q 


TV     3J77 

R     2.61  X  IO"   [f 


for  f  =  10  GHz 


Q,  *  14.  4x10 


Q,  is  sometimes  defined  as; 
f 


Q 


Af 


where  f  =  resonant  frequency  of  the  cavity 
Af  =  change  in  frequency  "between  half 
power  points 

Then  the  Af  for  this  mode  is  approximately  equal 


to  695  KHz 


2.3  Scattering  by  a  Spherical  Shell 

To  study  the  scattering  behavior  of  a  spherical 
shell,  it  will  be  necessary  to  expand  the  solution  of  the 
wave  equation  in  section  I  to  include  variations  with  0. 
This  expression  or  complete  solution  of  the  wave  equation 


13 


in  spherical  coordinates  is  given  in  terras  of  spherical 
vector  wave  functions  by  Stratton  (1941).   There  is  no 
loss  in  generality  to  assume  that  we  have  a  plane  wave 
incident  on  the  sphere  in  the  negative  z  direction  with 
the  electric  field  vector  in  the  positive  x  direction,  in 
a  spherical  coordinate  system  ( r,  6,      <t>)  related  to  the 
Cartesian  coordinate  system  by  x  =  r  sin  9   cos  <t>,  y  = 
r  sin6  sin  <t>,  and  z  =  r  cos  0.   Making  these  assumptions 
the  incident  field  equations  can  then  be  written  as: 


E  =  a   E   e     e 
x   o 


H  =  -a   En  e     e 

y  _£ 


jwt 
Normalizing  the  E  vector  to  unity  and  suppressing  the  e    , 

we  have : 


E  =  a   e 
x 


jkz 


H  =  -a   Y   e 
y  o 


jkz 


(3-1) 
(3-2) 


where  k  =  propagation  constant 

Z   =  intrinsic  impedance  of  free  space 

Y   =  intrinsic  admittance  of  free  space 
o 


14 
Employing  the  spherical  vector  wave  functions 
mentioned  above,  the  incident  field  is  expressed  in  spheri- 
cal coordinates  as: 

i    V^»   n  /  (2)       C1)  x      / 

E   =  X   j   2n  +  1       (Mnin  "  J  Neln)      (3"3) 
7  n(n  +  l)       01n       eln 

n=l 


00 


H1  =  jY   ^  /  2n  +  1  (IT    -  j  M(    )(3-4) 

0  7             n(n  +  1)  Oln  eln 

n=l 

(1)       (!) 

where  M     and  N     are  the  spherical  vector 

wave  functions  as  defined  by  Stratton  (1941)  as  given  below 


(!)    "    or  /   s   pm  ,      x   sin 
Me    =  +  m  ¥n(kr)   Fn  (c°s  q)  m<t>   a 

omn  kr         sin  $         cos 


m  ■>  s 

md>   a 


(3-5) 


-  *n(kr)    cK  (cos  gj   C0 
kr  de  sin 

(-1-)  .        ^m  cos 

N     =  n(n  +  1)   tn(kr)   Pn   (cos  0)      m<t>  a 
o"111  (kr)*  sin      r 

+  ^A(kr)   ^[pn   (cos  9)]     °°*    m  0>   a 

l  S  bill  D 

kr  de 


;   m  ^(kr)   Pn  (cos  0)   Sin  m«,   ^        (3-6) 
kr         sing      cos 
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where 


^n(kr) 


kr  J  (kr),  J   (kr)  is  the  spherical 
Bessel  function  of  order  n  and  the 
prime  indicates  differentiation  with 
respect  to  the  entire  argument. 


m 


P    (cos    0)    =    A    Legendre    polynomial 
n 

Substituting    equations    (3-5)    and    (3-6)     into 
equations    (3-3)    and    (3-4)    the    incident    fields    are    given 
as    follows: 


E      =   X       j       2n    +    1 
7  n(n    +    1 


tn(kr)         _gn       (cos    e)    cos     $  a 


kr 


sine 


n=l 


-*n(kr)        h\FJ    (cos    9)]        sin    <D 


kr 


Ze 


Pi 


-  j    n  ( n    +    1 )    V(kr)       Pn    (cos    9 )       cos    <t> 
(kr)2 

-j    ^n    (kr)  d[pn    (cos    9)]     c 


os    <p      a 


kr 


^9 


+j    ^h    (kr)       Pn    (cos    9 )       sin    0  a 


(kr) 


sin    9 


(3-7) 
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H   = 


n 

jy0  ST* J      2n  +  1 

/     n(n  +  1 


n(n  +  1)  ^n(kr)  Pn  (cos  q)     sin  <j> 
kr 


n=l 


+  tn(kr)     d[fn  (cos  0)]    sin 


$  a 


kr 


1 


^n(  kr )   Pn  (cos  B)       cos  <j>  a 
kr       sin  q 

+  j  ^n(  kr)   Pn  (cos  0)   sin  <J>  a 
kr       sin  0 

+  j  ^n(kr)   ^LPn(cos  flj    cos  0  * 
kr         de 


(3-8) 


Then  breaking  the  incident  fields  into  their 
respective  components,  we  have: 


E       =      -  j    cos     <j> 
(kr)2 


Y^    (2n    + 


>1 


1)       ^n(kr)  Pn(cos    B) 

kr 


(3-9) 


E 


00 

1 


n=l 


n 


cos    $  X     j       2n    +    1 
kr  n(n+l) 


n=l 


E      = 


00 

'I 


n 
sin    <t)\  "j       2n   +    1 

kr        /         n(n    + 1) 


*n(kr)    PnCcose)    -j    ^n(kr)    dfcj;    (cose)] 
sin    9  be 

(3-10) 

j    tn(kr)    Pjjfcos    0)-»n(kr)    o{  Pn(cosQ)] 
sin    g  ^0 

(3-11)" 


n=l 


H 


l   _ 


j'Y   sin  cf> 
o  r 


(kr) 


00 

I 

n=l 


j  ( 2n  +  1)   ^n(kr)   pJ(cos  9) 
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(3-12) 


H 


00 

1  v*  n       r.  t 

=  JY   sin  d>  \  J   2n  +  1   ^n 

kr    /  .    n(n+l)  _ 
n=l 


(kr)  d  [?n  (cos  9)]  + 


J  *n(kr)  Pn  (cos  9) 
sin  9 


(3-13) 


H   = 


I 


n 


j'Y   cos  <t>     y    j   2n  +  1 

kr   ^^    n(n  +  1) 

n=l 


"j  ^n(kr)  &[  P^(cos  9)] 


^9 


^(kr)  Pn  (cos  9) 
sin  9 


(3-14) 


The  scattered  field  can  be  represented  in  the  same 
manner  as  the  incident  field  with  spherical  vector  wave 
functions  as  follows: 


s 
E   = 


(4)  (4) 

U   M  _   -  j  B„  N  .   ) 
n   oln       n   eln 


(3-15) 


U) 


(4) 


H  =    ttQ)      (An  Noln  -  d  En  Meln 


(3-16) 


n=l 
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Where  M     and  IT  '    are  the  same  as  IVP  '    and  N^  '  except 


that  to  satisfy  the  requirement  that  the  scattered  field 

represent  an  outgoing  field  at  infinity,  ^n'kr)  is  replaced 

by  p  (kr)  =  kr  n\       (kr)  and  H  ^ c '    is  the  spherical  Hankel 
n  "  n 

function  of  the  second  kind  and  A   and  B   are  unknown 

n       n 

amplitude  coefficients. 

The  boundary  conditions  that  the  tangential  elec- 
tric fields  at  the  surface  of  the  sphere  must  equal  zero 
provide  a  means  for  solving  for  the  unknown  amplitude  coef- 
ficients A   and  B   as  follows: 

n       n 

T-.t    „i    „s 
E   =  E   +  E 


at  r  =  a 


E1  +  Es  =  0 


E1  +  Es  =  0 


then 


0  =  cos 


oo 

I 


n 


J   2n  +  1 


^^g    n(  n  +  1 

n=l 


♦n(ka)  Pn  (cos  e)    -   (3-17) 
sin  9 


-Pl 


3  *n  (ka)  a  [Pn  (cos  gj 

be 


+  cos  J^K  Pn(ka)  Pn  (cos  9)  -J  Bn  pn  (ka)  a[pl(cos  ej 
ka  /^  sin  9  $Q 


n=l 


j    in    (ka)    Pn    (cos    9)    - 
sin    0 
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(3-18) 


i     (ka)     d  [.  Pn (  c o s    6  )] 


+    sin 


d    B  p'     (ka)    Fn    (cos    9) 

n      n  ~ 

ka     Jm~*\  sin    0 
n=lL 


A       p    (ka)     5    [pn(  cos  9)] 

^9 


n      n 


Solving    equations    (3-17)    and    (3-18)    for    A      and 


B      we    get 
n 


An    =    -j       2n    +    1 


*J0    (ka) 


(3-19) 


n(n    +    1)  p-n    (ka) 


Bn    =    -jn    2n    +    1  tn(ka) 

n(n    +    1)  Pn    (ka) 


(3-20) 


The  scattered  transverse  electric  field  components 
for  large  kr  can  be  approximated  "by  replacing  p  (kr)  and 
p   '  (kr)  by  the  leading  terms  of  their  asymptotic  forms 
as  shown  by  Senior  and  Goodrich  (1964)  and  the  results  are: 
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„s    .         -jkr 
E   =  j  cos  4>  e 

e        ~^7 


I 

n=l 


j       2n  +  1 
n(n  +  1) 


*n  (ka)   S[Pn(cos  9)] 


LPn  (ka) 


*9 


(3-21) 


^n(ka)   Pn(cos  9) 


Pn(ka)     sin  9  _ 
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ES  = 


3    sin 

kr 


jkr 


I 


n=l 


jn  En  ±  1 


n(n  +  1) 


>n(ka) 


P  (ka) 


tn(ka)  ^[P^cos  9)1 
Pn(ka)      ^e 


P1(cos  0) 
sine 


(3-22) 


Er°   =  0 


A  more  concise  way  of  writing  the  equations  of  the  scattered 
fields  is: 


-Jkr    s 
j  cos  d)  e       s  (  8  ) 

kr  1 


(3-23) 


-j  sin 
kr 


-  jkr    ,, 

e       3S(e) 
2 


(3-24) 


s  s 

where:    S  (o)    and   3„(0)  are  the  summations  in 
1  2 

equations  (3-21)  and  (3-22). 


2.4  Scattering  by  a  Spherical  Shell  With 
a  Circular  Aperture 


This  problem  is  separate  from  the  one  in  the 
preceeding  section,  but  it  is  related  to  it.   As  we  have 
shown,  the  incident  and  scattered  fields  can  be  represented 
in  spherical  vector  wave  functions  as  follows: 


n  ,  (1)       (D, 

j   2n  +  1  (li  ■    N. 

— -, r-r     Oln       eln 

n(n  +  1) 


H   =  d 


E   = 


2 

n=l 

00 

iat    n(  n  + 
n=l 

00 

I 


1     (N     -  j  M  .  ) 
— —   v  Oln    °       ein; 
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(4-1) 

(4-2) 


(4)  (4) 

(An  M01n  "  3  Bn  Neln) 


n=l 


Hb  =  JY, 


00 

I 

n=l 


,(4 


(4=) 


An  N01n  ~  J  Bn  Meln 


(4-3) 


(4-4) 


And  applying  the  same  technique  the  fields  inside  the  shell 
can  be  represented  as: 


E 


00 

ln   =     >     [Cn  M01n    "    *    Dn  Neln] 


n=l 


H 


in 


n=l 


r  (H  (lV| 

[Cn  N03n  -     d    Dn   MelnJ 


(4-5) 


(4-6) 


The  far  field  scattering  coefficients  are  now:   Chanj 
and  Senior  (1969) 

00 

Ss     (0)    =    y      /[An    Pn    (cos    9)    -    Bn        a[Pn(cos    G)] 
•■?  sin    e  <$e 


S2    <9>    = 


00 

y       j  nUn    ^n    (cos    9)]  -    Bn       Pn(cos    6) 

4*mmti  be  sin    0 

n  =  l  L  J 


(4-7) 


(4-8) 
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And: 


-j  sin  ♦  exp( - jkr )   Ss   ( 
kr 


(4-9) 


Es   =  j  cos  ♦  exp( -jkr)    Ss  (0) 
0  J- 

kr 


(4-10) 


And  from  Senior  and  Goodrich  (1964) 


cr(e,  <t>)  =  ae(e,  *)  +  o-^e,*) 


(4-11) 


where 


o-Q(e,*) 


=   *•   cos   <t> 


Sls  (0) 


(4-12) 


a  (0,  <D)  =   *•    sin 


s2  <e> 


(4-13) 


Equations  (4-1)  through  (4-8)  show  that  once 

the  amplitude  coefficients  A  ,  B  ,  C  ,  and  D   are  deter- 

n'   n    n        n 

mined  the  fields  everywhere  are  known.  Therefore,  the 
problem  is  a  boundary  value  problem  where  the  boundary 
conditions  are  the  tangential  electric  fields  are  zero 
on  the  surface  of  the  shell  and  the  fields  through  the 
aperture  must  be  continuous.  When  these  boundary  condi- 
tions are  applied  to  the  equations,  it  is  possible  to 
solve  for  the  coefficients. 


23 
Chang  and  Senior  (1969)  discuss  a  number  of  ap- 
proaches to  this  problem  for  example,  using  integral  equa- 
tions in  terms  of  the  unknown  aperture  field,  and  the 
variational  technique  which  relates  the  backseat tering 
cross  section  to  a  scalar  function  U.   U  is  in  terms  of 
the  incident  electric  field  and  the  total  surface  current 
density  of  the  shell.   Both  of  these  approaches  require 
the  choice  of  a  set  of  expansion  functions^  for  the  integral 
equation  approach  a  set  of  expansion  functions  is  chosen  to 
represent  the  aperture  field  and  for  the  variational  method 
a  set  of  expansion  functions  is  chosen  to  represent  the 
surface  current  density.   The  accuracy  of  the  solution  for 
both  cases  is  determined  by  the  proper  choice  of  these 
expansion  functions.   Due  to  the  presence  of  the  cavity 
in  this  problem  the  proper  choice  would  be  extremely  dif- 
ficult. 

Due  to  the  difficulties  mentioned  above,  Chang 
and  Senior  (1969)  applied  the  method  of  the  least  square 
error  and  then  modified  this  method  to  separate  out  the 
discont inut ie s  at  the  edges  of  the  aperture.   This  ap- 
proach will  now  be  summarized:   (Refer  to  figure  2- 4a 
for  the  geometry  for  this  solution). 
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Applying  the  boundary  conditions  to  equations 

(4-1)  through  (4-6)  and  solving  for  the  coefficients  A 

and  B   in  terras  of  C   and  D   we  get: 
n  n       n 


A   p  (ka )  = 
n   n 


C   -  in  2n  +  1 
n       n(n  +  1) 


*n(ka) 


(4-14) 


Bn  pn(ka) 


D„  -  J 


2n  +  1 
n(n  +  1) 


\|f  (ka) 
n 


(4-15) 


Making  these  subsitut ions,  the  equations  for  the  scattering 


coefficients  A   and  B   are 
n       n 
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n=l 


for  9   <  9    <   n 
o     — 

And; 


\ a[pn  (cos  e?]-j    Bn    pn  (cos  ^ 
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for  0    <   6    <   9 
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where 


G1(e)  = 


00 

_\^  jn  2n  +  1 


1) 
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"        pl  (cos  e) 

i|r  (ka)  _2 

n 


sine 


-J  *B(ka) 


dpi  (cos  e)J 
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(4-2Q) 


(4-21) 


Applying  the  method  of  least  square  error,  Chang 
and  Senior  (1969)  obtained  a  set  of  2M  X  2M  matrix  equations 
which  when  solved  give  an  approximation  for  the  2M  scat- 
tering coefficients  A   and  B„  but  due  to  the  discontinuities 

n       n 

at  the  edge  of  the  aperture  these  equations  converge  very 
slowly.   To  improve  the  solution  Chang  and  Senior  (1969) 
examined  equations  (4-16)  through  (4-19)  and  separated  out 
the  discontinuities.   Then  applying  the  least  square  error 
method  to  the  new  equations  they  arrived  at  a  linear  system 
of  2(M  +  1)  equations  with  unknowns  defined  by: 


A 


n  =  an  +  Ki  *n(ka)  7n 


Bn  =  bn  +  K?  r= &n 

n     n     2  Pfl'(Ka)   n 


where 


Yn  -  Re  {cn} 


&n  =  Im  {cn} 
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c   = 


n    (  2n  -  l)(2n  +  3)   n    ( 2n  +  5 ) ( 2n  +  3 )   n+ 2 


( 2n  -  3)( 2n  +  1)   n" 2 


dn  =  exp  p(n  +  1/2)0O] 


which  converge  more  rapidly  giving  a  more  accurate  approxi- 
mation for  the  scattering  amplitude  coefficients  for  a 
given  M  than  the  original  least  square  method. 

The  results  of  this  "modified  least  square  method" 

produce  the  following  for  far  field  scattering  coefficients: 
M 
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(4-23) 


X 


2n 


=  -3 


**  (ka)   n 
yn 


X2M+1  ~  Kl 


X2M+2  "  'd  K2 


The  complete  development  of  this  approach  and 

the  computer  programs  to  carry  out  the  calculations  are 

presented  in  Chang  and  Senior  (1969).   The  theoretical 

response  for  an  aperture  size  of  0   =  30°  with  0=0, 

o 

0=0,  and  0.8  <  ka  <  4.8  with  M  =  10  for  ka  <  1.1,  M  =  20 
for  1.1  <  ka  <  4.8  and  M  =  30  in  the  region  of  expected 
peaks  as  calculated  by  Chang  and  Senior  (1969)  is  shown 
in  figure  2- 4b. 


29 


*-© 

"""— • 

o*> 

©— 

c 


J 


'•-^ 


.-©— 


^v 


•"»■» 


3 


.---'■ 


A' 


o 

Lf) 


o 


C 
CO 


o 

CM 


o 
o 


o 

00 


o 

c£> 


o 

CO 


<u 
oo 


o 

<u  -a 
-c  o 

Q..C 

oo  +-> 
a) 

4-  CD 
O    S_ 

rtJ 
CO    3 

c  cr 
o  CO 

+j   4-> 

o  to 

CD    fO 

oo  cd 

_J 

CO 

CO  "O 
O    CI) 

5-  •!- 

O  cj_ 

e  o 
•r-  s: 

J- 

CD     CD 

+->  -c 
-M  +-> 
<o 
o  ot 

CO    E 

.*:  i— 

CJ    CO 

to  rD 

DO 

"O  o 
CD  O 
+J  CO 

13 

Q-  II 
E 

O       O 
O   CD 


I 

CXI 

CD 

S- 

en 


,^/(0)° 


30 


2.5  Theoretical  Case  of  A  Remote  Cavity 


The  theory  in  the  preceeding  sections  will  now  be 
applied  to  a  remote  cavity.   The  assumption  is  made  that 
the  Q,  of  the  cavity  with  an  aperture  is  the  same  as  the 
closed  cavity.   It  is  known  that  the  aperture  in  the  cavity 
will  lower  the  Q,,  hut  the  calculations  employing  the  per- 
turbation theory  as  in  Slater  (1950)  are  extremely  compli- 
cated even  for  a  small  perturbation  and  are  beyond  the 
scope  of  this  thesis. 

The  fact  that  the  back  scattering  is  reduced  in 

the  vicinity  of  resonance  will  now  be  employed  to  determine 

the  detectable  change  in  frequency  and  subsequently  the 

change  in  N  units.   Assume  that  the  resonance  frequency  of 

the  cavity  is  10  GHz  (a  =  2.15  cm). 

For  a  change  of  100  N  units,  Af  =  f   -  10  GHz  =  999901 

'         0    1.0001 

Or  Af  for  a  change  of  one  N  unit  is  approximately 

10  KHz. 
Assuming  that  a  change  of  one  half  power  can  be  detected 
accurately  leads  to  the  assumption  that  a  change  of  one 
half  of  the  frequency  between  half  power  points  can  be 
detected.   Using  the  theoretical  maximum  Q,  for  the  TM 
mode  from  section  II  of  1.44  x  104  we  get  ££.   =    347.5  KHz 
and  represents  a  change  of  approximately  35  N  units. 


CHAPTER    III 


EXPERIMENTAL  RESULTS 


To  confirm  the  theory  in  the  preceding  sections, 
an  experimental  study  was  conducted. 

3.1   The  Model  and  The  Experimental  Set-up 

A  hollow  spherical  copper  shell  with  a  2  inch  out- 
side diameter  and  a  single  circular  aperture  of  60  =  15°, 
Figure  2-4a,  was  employed  as  the  resonant  cavity  model  to 
investigate  backscatter ing  near  resonance  from  such  a  cavity. 
Due  to  the  shape  of  the  field  lines  coupled  to  the  inside  of 

the  resonant  cavity  for  various  modes,  the  TE     mode  was 

J  '  101 

selected  because  the  circular  aperture  causes  less  distor- 
tion of  the  field  lines  for  this  mode.   The  theoretical  free 
space  (vacuum)  resonant  frequency  calculated  using  the 
outside  diameter  and  neglecting  the  aperture  was  8.45  GHz. 
(Found  later,  experimentally,  to  be  8.55  GHz.) 

An  anechoic  chamber  6  feet  wide,  8  feet  high  and  9 
feet  long  shown  in  Figure  3-la  was  constructed  to  eliminate 
extraneous  reflections  and  outside  interference.   The  resonant 

31 


32 


Figure  3-la.  Photograph  of  Experimental  Setup 
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spherical  cavity  model  to  be  examined  was  placed  on  a 
styrofoam  pedestal  inside  the  chamber  three  feet  from  the 
illuminating  antenna. 

The  equipment  arrangement  for  measuring  the 
backscattering  is  indicated  in  Figure  3-lb.   The  HP  624B 
signal  generator  was  used  as  the  transmitter  and  the  TS 
147D  signal  generator  was  employed  as  the  local  oscillator. 
By  frequency  modulating  the  TS  147D  with  the  sweep  genera- 
tor, the  requirement  for  continuous  frequency  tracking  of 
the  transmitter  was  eliminated.   Magic  tee  number  two 
served  as  both  the  mixer  and  detector.   The  1216A  ampli- 
fier was  used  as  an  IF  amplifier  and  the  output  read  di- 
rectly off  the  oscilloscope  face.   The  oscilloscope  response 
represented  the  band-pass  characteristics  of  the  IF  ampli- 
fier, but  the  peak  response  was  directly  related  to  the 
amplitude  of  the  electric  field  reflected  from  the  spher- 
ical shell. 

The  procedure  for  taking  a  measurement  was  as 
follows:   With  the  transmitter  tuned  to  the  desired  fre- 
quency and  the  spherical  cavity  model  removed  from  the 
anechoic  chamber,  arms  A  and  B  of  the  first  magic  tee  were 
balanced  by  tuning  the  three  tuning  screws  in  arm  A  and 
the  variable  attenuator  in  arm  B  until  no  signal  was 
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observed  on  the  oscilloscope.   This  meant  that  none  of 
the  transmitter  signal  from  arm  C  was  entering  arm  D  and 
that  the  transmitter  signal  was  divided  equally  between 
arms  A  and  B.   The  spherical  cavity  model  was  then  placed 
on  its  pedestal  and  irridiated.   One  half  of  the  signal  re- 
flected from  the  model  entered  arm  D  and  was  mixed  and 
detected  in  magic  tee  number  two.   The  detected  signal  was 
then  amplified  and  recorded.   The  model  was  rotated  from 
its  initial  position  of  6  =  0°  to  6  =  90°  in  10°  steps  and 
the  measurements  repeated  for  each  step.   The  final  step 
was  to  rotate  the  model  to  6  =  180°  (circular  aperture  on 
the  back)  and  make  a  measurement, 

3.2   Measurements 


The  initial  measurements  were  made  to  determine 
the  resonant  frequency  of  the  cavity.   The  copper  shell 
has  a  finite  thickness  which  causes  the  cavity  to  be 
smaller  than  the  outside  measurement  would  indicate  and 
increases  the  resonant  frequency  while  the  aperture  has  a 
tendency  to  lower  the  resonant  frequency.   The  theory  shows 
that  the  signal  backscattered  from  the  aperture  (6  =  0°)  is 
less  than  from  a  solid  shell  (9  =  180°)  only  in  the 


36 


vicinity  of  resonance.   For  the  above  reasons,  measurements 
of  the  backscatter  from  the  aperture  and  the  solid  shell 
were  taken  starting  from  a  frequency  below  the  theoretical 
resonant  frequency  and  recorded  as  frequency  was  increased. 
The  experimental  resonant  frequency  was  found  to  be  approxi- 
mately 8.55GHz. 

After  the  resonant  frequency  was  determined, 
measurements  of  backscattering  were  made  in  the  vicinity 
of  this  frequency.   The  ratio  of  the  backscattered  signal 
from  the  shell  with  9=0°  to  the  backscattered  signal  from 
the  shell  with  ©  =  180°  was  recorded  and  plotted.   The 
first  attempts  to  obtain  repeatable  data  failed  and  was 
believed  to  be  due  to  the  inaccuracy  of  the  frequency 
measuring  device  on  the  transmitting  equipment  and  a 
considerable  amount  of  copper  oxide  present  on  the  copper 
shell  surface.   The  shell  was  cleaned  and  the  measurements 
repeated.   It  remained  extremely  difficult  to  obtain  re- 
peatable data,  but  an  attempt  was  made  and  the  results  of 
three  sets  of  data  were  averaged  to  obtain  the  curve  shown 
in  Figure  3-2a.   The  curve  in  Figure  3-2a  was  obtained  by 
computing  the  ratio  of  the  signal  backscattered  from  the 
shell  with  9=0°  to  the  signal  backscattered  from  the 
shell  with  9  =  180°  and  the  results  plotted  versus  frequency. 
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This  curve  provides  a  means  of  determining  the  approximate 
bandwidth,  the  resonant  frequency,  and  the  detectable 
change  in  N  units. 

As  the  backscattering  measurements  above  were 
being  made,  the  variations  in  the  backscattering  as  a  func- 
tion of  G  were  also  recorded  in  10°  steps  from  9=0°  to 
9  =  90° .   The  purpose  of  this  measurement  was  to  determine 
the  effects  of  the  aperture  as  a  function  of  illumination 
angle  much  as  what  would  happen  in  the  use  of  a  free  fall- 
ing remote  cavity.   The  results  of  these  measurements  are 
shown  in  Figures  3-2b  through  3-2d.   These  curves  were  ob- 
tained by  taking  the  ratio  of  the  signal  backscattered  from 
the  shell  at  an  angle  9  to  the  signal  backscattered  from 
the  shell  at  9  =  180°  with  the  frequency  held  constant  and 
the  results  plotted  as  a  function  of  9.   Figures  3-2b 
through  3-2d  show  the  variations  in  the  backscattering  as 
a  function  of  9  for  three  different  frequencies,  Figure 
3-2b  was  plotted  for  a  frequency  slightly  below  resonance 
while  the  data  for  Figure  3-2d  was  taken  at  a  frequency 
slightly  above  resonance  and  the  data  for  Figure  3-2c  was 
taken  at  the  approximate  resonant  frequency  of  the 
c  avity . 
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3.3   Discussion  of  Experimental  Results 

The  objectives  of  the  experiment  were  to: 

(1)  Determine  the  resonant  frequency  of  the  cavity 
with  a  circular  aperture. 

(2)  Determine  the  band  width  of  the  cavity  with  a 
single  aperture. 

(3)  Determine  the  change  in  N  units,  of  the  gas  fil- 
ling the  cavity,  that  could  be  reliably  detected. 

(4)  Determine  the  effects  of  angle  of  incidence  on 
the  copper  sphere  relative  to  the  location  of 
circular  aperture. 


The  resonant  frequency  was  found  to  be  approxi- 
mately 8.55  GHz  and  the  backscattering  was  investigated  in 
the  vicinity  of  this  resonant  frequency.   Figure  3-2a  re- 
presents the  results  of  this  investigation.   From  the  curve 
in  Figure  3-2a,  the  approximate  "Q"  of  the  cavity  with  a 
single  aperture  was  found  to  be  1600  with  a  corresponding 
Af  of  approximately  5.4  MHz.   This  Af  would  allow  the  de- 
tection of  a  minimum  change  of  approximately  320  N  units. 
Since  the  maximum  change  in  N  units  in  the  atmosphere  from 
a  vacuum  to  sea  level  is  approximately  340,  this  system 
would  reveal  essentially  no  useful  information  about  the 
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index  of  refraction  of  the  atmosphere  as  a  function  of 
height.   To  detect  a  change  of  one  N  unit,  the  "Q"  of  the 
cavity  would  have  to  be  505,000  or  have  a  Af  of  approxi- 
mately 17KHz  at  a  resonant  frequency  of  8.55  GHz. 

The  most  difficult  part  of  the  experiment  was 
that  of  obtaining  repeatable  data.   This  was  due  in  a 
large  part  to  the  inaccuracy  of  the  frequency  measuring 
device  on  the  HP  624B.   The  accuracy  of  the  frequency 
measuring  device  was  ±  0.03  percent  or  approximately 
_+2.565MHz  and  this  error  could  be  as  much  as  one  half  Af, 
the  bandwidth.   This  represents  a  significant  error,  how- 
ever, the  frequency  measurements  were  accurate  enough  to 
definitely  show  that  the  bandwidth  of  the  resonant  cavity 
was  much  greater  than  the  bandwidth  needed  to  be  able  to 
detect  a  change  of  one  N  unit  in  the  dielectric  filling  the 
cavity.   Another  error  was  introduced  in  the  measurements 
by  a  noise  level  in  the  measuring  system  of  0.4V.   This 
created  a  problem  only  when  the  desired  signal  level 
dropped  to  this  value  and  really  had  very  little  effect 
on  the  measurements. 

To  satisfy  the  fourth  objective  of  the  experiment, 
the  measurements  for  Figures  3-2b  through  3-2d  were  taken. 
Investigating  Figure  3-2c,  we  see  that  there  is  little 
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effect  caused  by  a  change  in  the  illumination  angle  9  of 

±  20  degrees  relative  to  the  location  of  the  circular  aper- 
ture but  a  considerable  effect  for  changes  in  9  of  more 
than  20  degrees.   Figure  3-2b  shows  that  the  changes  of 
the  backscatter ing  with  changes  of  9  are  essentially  the 
same  when  the  frequency  of  the  probing  signal  is  slightly 
below  resonance  as  when  it  is  at  resonance.   However,  when 
the  frequency  is  slightly  above  resonance  the  changes  are 
very  different  from  those  at  resonance. 

Measurement  errors  that  could  be  introduced  by 
inaccurate  positioning  of  the  model  are  considerable,  for 
example,  if  the  spherical  copper  shell  had  been  rotated 
to  an  angle  of  0  ±  25°,  Figure  3-2d,  the  resonant  fre- 
quency would  appear  to  be  2  MHz  higher  than  the  actual 
resonant  frequency  which  would  represent  a  change  of  234  N 
units  in  the  dielectric  filling  the  cavity.   On  the  other 
hand,  if  the  shell  had  rotated  more  than  ±  60°  there  would 
be  no  change  in  the  backscatter ing  for  small  changes  in 
frequency . 


CHAPTER  IV 


SUMMARY 


The  feasibility  of  remote  measuring  the  index  of 
refraction  of  the  atmosphere  using  the  resonance  charac- 
teristics of  a  hollow  spherical  cavity  with  apertures  was 
investigated  both  theoretically  and  experimentally.   The 
theoretical  study  assumed  a  closed  cavity  and  produced  a  Q, 
with  a  bandwidth  narrow  enough  to  enable  the  detection  of 
a  change  in  the  index  of  refraction  of  the  gas  filling  the 
cavity  of  approximately  35  N  units.   The  experimental 
study  produced  a  minimum  detectable  change  of  approximately 
320  N  units. 

To  produce  a  feasible  system  for  measuring  the 

index  of  refraction  of  the  atmosphere  a  minimum  detectable 

change  of  no  more  than  1  N  unit  would  be  required.   As 

shown  by  both  the  experimental  and  theoretical  studies, 

the  sensitivity  of  the  proposed  system  using  a  hollow 

spherical  resonant  cavity  with  one  aperture  is  far  from 

that  required.   In  actual  practice  a  number  of  apertures 

would  be  required  for  adequate  flushing  of  the  cavity  and 

to  make  the  measurements  less  dependent  on  the  angle  of 

illumination,  Q. 
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With  more  accurate  transmitting  equipment,  the 
measurements  could  be  repeated  and  plotted  more  accurately, 
but  since  the  equipment  used  in  this  experiment  show  con- 
clusively that  the  bandwidth  of  the  cavity  is  far  too 
large  to  enable  the  measurement  of  a  minimum  value  of  N 
with  only  one  aperture,  it  is  not  recommended  that  this 
method  of  measuring  the  index  of  refraction  be  further 
explored . 
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